Molecular gas in extreme star-forming environments: the 
starbursts Arp 220 and NGC 6240 as case studies 
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ABSTRACT 

We report single-dish multi-transition measurements of the 12 C0, HCN, and 
HC0 + molecular line emission as well as HNC J=l-0 and HNCO in the two ul- 
traluminous infra-red galaxies Arp 220 and NGC 6240. Using this new molecular 
line inventory, in conjunction with existing data in the literature, we compiled the 
most extensive molecular line data sets to date for such galaxies. The many ro- 
tational transitions, with their different excitation requirements, allow the study 
of the molecular gas over a wide range of different densities and temperatures 
with significant redundancy, and thus allow good constraints on the properties 
of the dense gas in these two systems. The mass (~ (1 — 2) x 10 10 M ) of dense 
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gas (> 10 5 ~ 6 cm -3 ) found accounts for the bulk of their molecular gas mass, and 
is consistent with most of their IR luminosities powered by intense star bursts 
while self-regulated by 0,B star cluster radiative pressure onto the star-forming 
dense molecular gas. The highly excited HCN transitions trace a gas phase 
~ (10 — 100) x denser than that of the sub-thermally excited HCO + lines (for 
both galaxies). These two phases are consistent with an underlying density-size 
power law found for Galactic GMCs (but with a steeper exponent), with HCN 
lines tracing denser and more compact regions than HCO + . Whether this is 
true in IR-luminous, star forming galaxies in general remains to be seen, and 
underlines the need for observations of molecular transitions with high critical 
densities for a sample of bright (U)LIRGs in the local Universe - a task for which 
the HI-FI instrument on board Herschel is ideally suited to do. 

Subject headings: galaxies: starbursts - galaxies: individual (Arp 220 and NGC 6240) 
- ISM: molecules 



1. Introduction 

A cardinal result in extra-galactic astronomy over the last 20 years has been the discovery 
by IRAS of a population of local Ultra Luminous Infra- Red Galaxies (ULIRGs - Soifer et al. 
(1986); see also Sanders & Mirabel (1996), and Lonsdale, Farrah & Smith (2006)). Although 
these galaxies are rare in the local Universe, their extreme IR-luminosities (L m ~ 10 12 L & ) 
and merger morphologies strongly suggest that they signify a highly active and important 
phase in galaxy evolution - a scenario which is corroborated by the rapid increase in their 
space density with redshift. Furthermore, ULIRGs - being gas-rich mergers harboring the 
most extreme star-forming conditions encountered in the local Universe - are thought to 
re-enact the galaxy formation processes we are only barely able to glimpse in the distant 
Universe. Although, detailed studies of the interstellar medium (ISM) in proto-galaxies 
at high redshifts will be possible with the new generation of cm/mm/sub-mm telescopes 
and receivers coming online within the next decade, our expectations of how such objects 
will look in mm/sub- mm line emission are ill-informed due to the dearth of similar data 
locally. Multi-line studies of the low-intensity starbursts M82 and NGC 253 have been made 
( Nguyen- Q-Rieu et al. 1992; Jackson et al. 1995; Paglione, Jackson & Ishizuki 1997), but 
given that their star formation rates are about three orders of magnitude lower than that 
of typical ULIRGs, they may very well not be good guides for the state of the ISM in local 
ULIRGs and high-z starbursts (Aalto et al. 1995; Papadopoulos & Seaquist 1998; Downes 
& Solomon 1998). 
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In contrast, local ULIRGs, owing to their extreme IR/CO-brightness, are ideal targets 
for multi-line studies of heavy rotor molecules such as HCN, CS and HCO + . These lines 
are usually 5 — 30 times weaker than their CO J-level counterparts, making them currently 
impossible to detect in more distant galaxies and/or galaxies with lower star formation rates. 
Secondly, their compact H2 distributions - as revealed by interferometric maps (e.g. Scoville 
et al. 1997; Tacconi et al. 1999) - allow us to collect the total line flux of these often weak lines 
with a single pointing. This is important since multi-beam receivers, capable of observing 
many transitions of a variety of species within a reasonable amount of telescope time, are still 
lacking. These two factors, brightness and relatively compact molecular gas distributions, 
allow us to study the properties of the dense (n(H 2 ) > 10 4 cm -3 ) gas phase in local ULIRGs in 
much greater detail than would otherwise be possible. This is of great importance, because 
the dense gas provides the immediate star formation fuel in galaxies (Gao & Solomon 2004a,b; 
Wu et al. 2005). The possibility of a universal star formation efficiency, as proposed by these 
studies, can now be checked in a much more powerful manner, since we can estimate the 
dense H2 mass using multiple lines rather than a single line such as HCN J = 1 — 0. 

In the local Universe, Arp 220 and NGC 6240 are ideally suited for studies of the dense 
gas in extreme starburst galaxies, as they are the two most studied (U)LIRGs to date. 
They have the best sampled FIR/sub-mm spectral energy distributions (SEDs) of their class 
(e.g. Dopita et al. 2005), and the largest number of mm/sub-mm CO, HCN, CS, and HCO + 
transitions already observed. Thus not only are their dust continuum SEDs well constrained, 
but it is also relatively straightforward to fill the gaps in their molecular line inventory. Being 
the nearest and best studied, Arp 220 is the ULIRG most frequently used as a template for 
high-z studies of starbursts. The prominent Active Galactic Nucleus (AGN) present in 
NGC 6240 (even double- AGN, Komossa et al. (2003), but not in Arp 220, cf. Iwasawa et 
al. (2005)) offers the possibility to study the influence of an AGN on the molecular ISM, 
for example how X-ray Dissociation Regions (XDRs) affect H 2 cloud chemistry, significantly 
altering the HCO + and HCN abundances (Lepp & Dalgarno 1996; Usero et al. 2004). 

In this paper we present single-dish observations of a number of molecular emission lines 
from Arp 220 and NGC 6240, which together with results from the literature allows us to 
compile a very large (and almost identical) molecular transition catalog for these two extreme 
starbursts. Armed with this unique catalog, and the large range of excitation conditions it 
can probe (with significant overlap/redudancy between HCN and CS lines) we use it to: 
a) set constraints on the physical properties of the dense molecular gas phase, b) find the 
gas density probed per molecular species and its mass range, and c) explore underlying star 
formation efficiencies, density-size hierarchies, and examine the species best suited for tracing 
the true star forming phase. 
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The paper is outlined as follows: observations and results are described in £J2] and [31 
respectively. §H details the analysis of our observations, which includes deriving dense gas 
mass and star formation efficiencies based on the radiative transfer modeling of the observed 
lines. In §[5] we discuss the implications of our findings in the context of upcoming missions 
such as Herschel. In the flat cosmology (O m = 0.27, Q\ = 0.73, and H = 71 kms -1 Mpc -3 - 
Spergel et al. 2003) adopted throughout this paper the luminosity distances to Arp 220 and 
NGC6240 are 77.5 and 105.5 Mpc, respectively. 



2. Observations 

The 15-m James Clerk Maxwell Telescop^j] (JCMT) on Mauna Kea, Hawaii was used 
to observe a number of 12 CO, 13 CO, HCN, CS, and HCO+ lines in Arp 220 and NGC6240 
(Table [[]). These observations were done in good (t225gh z < 0.035) to medium (0.08 < 
r 225GHz ^ 0.12) weather conditions during several runs on the JCMT between January 2003 
and December 2004. The two JCMT receivers A3 (2 11-279 GHz) and 53(312-370 GHz) were 
employed in order to observe all the lines, with the B3 receiver tuned in SSB mode. In order 
to ensure baselines as flat as possible, the observations were done in fast beam-switching 
mode with a chop throw of 30" (AZ) and a chopping frequency of 1 Hz. As backend we used 
the Digital Autocorrelation Spectrometer (DAS) in its widest, ~ 1.8 GHz (~ 1100 kms" 1 
for the B3 receiver) mode to cover the full velocity-width of the expected broad lines. The 
pointing was checked every 1-2 hr throughout each observing run and was ~ 3'/5 (rms) or 
better. Aperture and main beam efficiencies were derived from observations of Mars and 
Uranus throughout the observing runs and were compared to mean JCMT values over for 
the observing periods. There is excellent agreement between the two sets of values so we 
adopt the latter under the assumption the mean values carry less inherent uncertainty than 
single measurements. The details of the observations, such as beam sizes, line frequencies, 
system temperatures and total on-source integration times, are summarized in Table [H 

The HCN 1-0 and 13 CO 2-1 lines were observed with the IRAM 30-m telescope on Pico 
Veleta (Spain) in June 2006 as part of a larger molecular line survey of LIRGs (Papadopoulos 
et al. 2007a). Each line was observed simultaneously in two polarizations using a pair of SIS 
receivers. The 3 mm A100/B100 receivers were used for the HCN 1-0 line in conjunction 
with two 1 MHz filter banks, each with 512 channels and a channel separation of ~ 3.4 kms -1 



1 The JCMT is operated by the Joint Astronomy Centre on behalf of the United Kingdom Particle Physics 
and Astronomy Research Council (PPARC), the Netherlands Organisation for Scientific Research, and the 
National Research Council of Canada. 
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at the observed frequency of 87 GHz. For the 13 CO 2-1 line, we used the 1mm A230/B230 
receivers together with the two wide 4 MHz filter banks, each with 1 GHz bandpass and a 
channel separation of ~ 5.5 km s -1 at the observed frequency of 216.5 GHz. The beam size 
of the 30-m telescope is 28"at 115 GHz and 12"at 217 GHz. The data were acquired with the 
New Control System (NCS) in series of four minute scans, each comprised of eight 30 s sub- 
scans. Total integration times were 1.2 hours for Arp220 and almost 2 hours for NGC6240 
in HCN 1-0 and 13 CO 2-1 (Tabled]). To obtain the flattest possible baselines, the wobbler 
switching (nutating subreflector) observing mode was employed with a frequency of 0.5 Hz 
and beam throws of 180 — 240". Both targets were observed under good weather conditions 
in the late night hours with T sys ~ 95 — 125 K at 3 mm and 225 — 500 K at 1 mm on the 
antenna temperature (TjJ) scale. Pointing and focus were checked frequently throughout 
observations, and the pointing errors were found to be ~ 3" (rms). 

Finally, we present previously unpublished observations of Arp220 obtained with the 
IRAM 30-m during the period 1988-1992. These data include 13 CO 1-0, CS 2-1, 3-2, C 18 
1-0, and HNC 1-0, the isomer of HCN. For these observations, SIS receivers were used 
together with the 512 x 1 MHz spectrometers. The observations were done with a wobbler 
switching of 180-240" in azimuth. The data were calibrated with cold and ambient loads 
and, depending on frequency and weather, single sideband temperatures were 200-600 K. 
The pointing, checked with planets and quasars on a regular basis, was ~ 3" (rms). 

The spectra of all the molecular lines observed in Arp 220 and NGC 6240 as part of this 
study are shown in Fig. [1] and [21 respectively. The spectra have been plotted on the T\ 
temperature scale, i.e., the antenna temperature corrected for atmospheric absorption and 
rearward spillover (Kutner & Ulich 1981). The velocity-integrated line flux densities were 
estimated using 



where r] ap is the aperture efficiency, D the telescope diameter, and ks Boltzmann's con- 
stant (Gordon, Baars & Cocke 1992). We find T = 15.625 JyK" 1 for the JCMT and 
T = 3.907 JyK -1 for the IRAM 30-m Telescope. T\ is the baseline subtracted spectrum, 
and K c (x) = x 2 /(l — e~ x2 ),x = 6> S /(1.2# HPBW ) (where s =source diameter) accounts for the 
coupling of the gaussian part of the beam with a finite-sized, disk-like source. Interferometric 
observations have shown that the CO emitting regions in Arp 220 and NGC 6240 are ~ 4" 
in angular diameter (Tacconi et al. 1999; Scoville et al. 1997). Our line observations collect, 
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therefore, all the flux, with minimal dependence on K c , which is at most ~ 1.04 for the 
highest frequency (and thus smallest beam) observations reported here. 

All of the measured line fluxes along with their errors are tabulated in Table [3J The 
errors include thermal and calibration errors, as well as errors on the assumed aperture 
efficiencies. Each of these are assumed to be uncorrelated and thus added in quadrature. 
The latter two sources of error were found to contribute by no more than 10 per-cent each, 
while the thermal noise was calculated for each line using: 

<*{ [ n dv) = JWZ (l + ^) 1 a(T* AtCh )Av ch , (3) 

where N& v (= Av/Av c h) is the number of channels across the line, A& as is the number of 
channels (each Av c h wide) used to estimate the baseline level, and cr(T\ ch ) is the channel- 
to-channel rms noise, assumed to be uniform across the spectrum. 

2.1. Calibration issues with the A3 receiver 

At v > 230 GHz the A3 receiver on the JCMT may present a calibration problem that 
has to do with the fact that the sideband ratio in this Double Sideband (DSB) receiver is 
not unity but can vary as a function of the LO frequency. This may have affected our HCN 
3-2 measurements, which we now reduce separately for each sideband, and intend to correct 
properly before co-adding. 

In the standard calibration the raw normalized spectrum is scaled by T ca \ = (l+Gj s )AT, 
where Gi S = Gi/G s (= 1 for DSB mode) is the ratio of the image and signal sideband gains, 
and AT depends on the mean atmospheric and telescope cabin temperatures, the telescope 
transmission efficiency, and the line-of-sight optical depth (see Eq. A12 of Kutner & Ulich 
1981). The first step for reducing spectra for which G,i s =/ 1 is to divide the spectra by 2 
in order to undo the standard online calibration. Next, the LSB and USB tuned spectra 
were averaged and corrected separately using the G,i s {vlo) dependence given on the JCMT 
website^], which for for the HCN 3-2 line yielded Gi S = 0.81. Finally, the resulting corrected 
spectra were co-added. 



' http://www.jach.hawaii.edu/JCMT/spectralJine/InstrumentJiomes/RxA3i/sideba 
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3. Results 

In this section we describe the new transitions observed in Arp 220 and NGC 6240 as part 
of this study, along with previously detected transitions that we have re-observed. In order 
to compare the line strengths not only between the two galaxies but also between different 
transitions within the same system, we have calculated velocity-integrated line luminosities 
for all of the observed transitions. These were calculated using 

U = 3.25 x 10 7 ( ^)~ 2 ( (1 + ( ^ S " dV 1 ) , (4) 

VGHz/ V M P C / VJykms- 1 / V; 

(Solomon, Downes & Radford 1992b). The resulting line luminosities calculated in this 
manner are tabulated in Table [3j 



3.1. Arp 220 

In Arp 220, high-resolution interferometric observations have revealed the double-horn 
profiles of the 12 CO 1-0 and 2-1 lines (Scoville et al. 1991) to originate from two nuclei 
embedded within a circumnuclear ring or disk (Scoville et al. 1997; Downes & Solomon 
1998). We see the same line profile in our 12 CO 2-1 spectrum (Fig. [Ik,), with a strong 
emission peak at ~ —200 km s" 1 (relative to Vlsr — 5454 km s -1 ) and a less intense peak 
at ~ +50kms~ 1 . These two peaks correspond, respectively, to the western, blueshifted 
nucleus and to the eastern, redshifted nucleus in the molecular disk (Downes & Solomon 
1998). Our 12 CO 3-2 line, single-dish spectrum show a double-horn profile similar to the 
lower transitions (Lisenfeld et al. 1996; Mauersberger et al. 1999). 

One subtle difference exists however, and is confirmed by our own 12 CO 3-2 spectrum 
(Fig. [Tb). The relative strength of the peak at ~ +50 km s _1 is greater in the 3-2 line than in 
the 1-0 and 2-1 spectra. In the 3-2 spectrum the western component is only slightly stronger 
than the eastern. Recent interferometric observations with the Submillimeter Array (SMA) 
showed that the CO 3-2 and the 890-/im dust emission is indeed strongest in the western 
nucleus (Sakamoto et al. 2008). 

The 13 CO lines observed in Arp 220 with the IRAM 30-m and the JCMT are shown in 
Fig. [lb and d, respectively. Although the 13 CO lines seem somewhat narrower than their 
12 CO counterparts, this may be simply an effect of the low signal-to-noise ratio in the line 
wings and the uncertain baseline level. In general all the 12 CO and 13 CO lines seem to have 
very similar line profiles, within the available signal-to-noise. Both the 13 CO 1-0 and 2-1 
spectra observed with the IRAM 30-m show tentative evidence of double peaks, consistent 
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with the 12 C0 lines. In the 1-0 spectrum, both peaks are equally strong, whereas the peak 
at ~ —200 km s" 1 is stronger in the 2-1 spectrum. This also appears to be the case in the 2-1 
and 3-2 spectra observed with the JCMT (Fig. [TH), although we note the 3-2 spectrum in 
particular is very noisy. The fact the 13 CO 2-1 spectra obtained with the two telescopes are 
in agreement in terms of line shape as well as flux density gives us confidence our observations 
are able to robustly discern the line shapes of even faint lines like 13 CO 2—1. 

The HCN 1-0, 3-2, and 4-3 spectra observed in Arp220 are shown in Fig. dfe-g, re- 
spectively. The HCN 1-0 spectrum shows a clear double-horn profile similar to 12 CO 3- 
2, with equally strong peaks. The HCN 3-2 and 4-3 line profiles are both very broad 
(fwhm ~ 540 — 590kms _1 ) and share the same overall shape. They both show evidence of 
double peaks with the ~ -l-SOkms" 1 component being more prominent. The shape is unlike 
the HCN(l-O) spectrum - as well as the 12 CO 2-1 and 3-2 spectra - where the ~ —200 km s _1 
component is brighter. 

The CS lines - of which only CS 3-2 has been previously reported for Arp 220 - trace 
a much overlaping excitation range with that of HCN (see Table 6), and are shown in Fig. 
[Th-k. All CS lines have double-horn profiles, with the ~ +50kms _1 peak stronger in the 
CS 3-2, 5-4, and 7-6 spectra, while the two peaks are almost the same strength in the 
2-1 spectrum. This mimicks the observed evolution from low- to high-J in the HCN line 
profiles. We note, however, the apparent asymmetry of the CS 3-2 line may be an artifact 
of the limited receiver bandwidth. In determining the baseline, we forced the right edge of 
this spectrum to zero, so our line flux estimate for this line is likely an underestimate. 

Finally, in the last column of Fig. [TJ we show the HCO + 4-3, HNC 1-0, and the blended 
C 18 1-0 and HNCO 5 ,5-4 ,4 lines. Similarly to the HCN and CS lines, the HCO+ spectrum 
peaks at ~ +50kms _1 , and shows some faint emission at ~ — 200 km s -1 , which corresponds 
to the stronger 12 CO feature. In this respect, the HCO + spectrum looks similar to CS 
7-6, where the western component is almost completely quenched. In stark contrast to 
this, the double-horn spectrum of HNC 1-0 shows a remarkably strong emission peak at 
~ — 200 km s^ 1 , suggesting a substantial fraction of the HNC is concentrated in the western 
nucleus. The blended C 18 1-0 and HNCO 5o,5~4o,4 spectrum is shown in Fig. [Th. In order to 
determine separate parameters for C 18 1-0 and HNCO 5o,5-4o,4, two gaussian components 
with the same width and fixed separation were fitted. From the resulting intensities we 
calculated integrated line fluxes and luminosities (Table [3]). 

It is telling that the double-horn structure seen in the low-J 12 CO spectra is also dis- 
cernible in the high-J spectra of HCN, CS, and HCO + , albeit with the relative strengths of 
the two peaks reversed. In fact, a rather clear trend is apparent: the ~ — 200kms _1 peak is 
stronger in the 12 CO 2-1 spectrum, the peaks are about equally strong in the 12 CO 3-2 and 
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HCN 1-0 profiles, and the ~ +50kms~ 1 peak is stronger in the HCN 3-2 and 4-3 spectra. 
Thus, as we climb the J-ladder and probe higher and higher critical densities, the strength 
of the emission shifts from the western to the eastern nucleus. This picture is confirmed by 
the CS 7-6, and HCO + 4-3 spectra which have even higher critical densities than the HCN 
1-0 line and are completely dominated by the ~ +50kms~ 1 component with very little or 
no emission at all at ~ — 200 km s -1 . 

An important point to make here is that it is only due to the high-density gas tracers 
such as HCO + , HCN and CS that we are able to conclude that the eastern nucleus harbours 
most of the dense gas. Without those, we would have been inclined to conclude, based on 
the CO 3-2 spectrum (Fig. [lb) and its distribution (Sakamoto et al. 2008), that most of the 
dense gas was contained in the western component. 

Given the modest (yet real) differences of the relative line strengths across the velocity 
range of Arp 220, and for reasons of simplicity and straightforward comparison to NGC 6240, 
we decided to use the line intensities average over its entire profile and the associated line 
ratios in our susbesquent analysis (§3J. 

3.2. NGC 6240 

The 12 CO 2-1 and 3-2 spectra of NGC 6240 are shown in Fig. [2^ and b, respectively. 
Unlike Arp 220, the 12 CO 3-2 and 2-1 spectra in NGC 6240 have almost identical shapes, 
with spectra having a single-peaked profile, skewed towards the blue (~ — 200kms _1 relative 
to vlsr — 7359 kms -1 ). This is consistent with high-resolution interferometric observations 
of low- J 12 CO lines in the inner regions of NGC 6240, which suggest more than half of 
the molecular gas resides in a central structure dominated by turbulent flows, rather than 
systematic orbital motion (Tacconi et al. 1999). The 13 CO 2-1 spectra obtained with the 
JCMT and IRAM 30-m are shown in Fig. [2fc and d, respectively, and both appear similar in 
shape to their 12 CO counterparts. Unfortunately, we cannot compare the 13 CO 3-2 spectrum 
(Fig. [2U) as no obvious detection of this transition was made. We are able, however, to set 
a sensitive upper limit on its integrated flux (Table [3]). 

The HCN 1-0, 3-2, and 4-3 transitions observed in NGC 6240 are shown in Fig. [2k- 
g. Although the HCN 1-0 line has previously been detected (Solomon, Downes & Radford 
1992a; Tacconi et al. 1999; Nakanishi et al. 2005), this is the first time the 3-2 and 4-3 
transitions have been detected in this galaxy. All three HCN lines peak at ~ —200 km s -1 
and are consistent with the CO line profiles. 

Unlike Arp 220, where several CS lines were observed, only the 7-6 line was observed 
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in NGC 6240 due to limited observing time. The spectrum is shown in Fig. [2h. This is the 
only detection of CS in NGC 6240 to date and the line profile is consistent with those of CO 
and HCN. 

Finally, the HCO+ 4-3 spectrum observed in NGC 6240 is shown in Fig. H. The HCO+ 
4-3 line profile is in good agreement with the 1-0 line (Nakanishi et al. 2005) and is consistent 
with a single peak at ~ — 200kms _1 . 

In summary, all of the transitions detected in NGC 6240 exhibit a similar, Gaussian- 
like velocity profile, peaking at about the same velocity (~ — 200 km s -1 ), in contrast to 
Arp 220 where all the lines not only show a double-peaked profile but a strong evolution in 
the dominant peak is also discernible. The similar velocity profile over such a large range 
of densities probed in NGC 6240 signifies a kinematical state of the molecular gas where the 
various phases remain well-mixed down to much smaller scales than in Arp 220 where the 
denser phases seem settled to two counter-rotating dense gas disks, enveloped by a much 
more diffuse CO-bright gas phase. Sub-arcsec interferometric imaging of the HCN, CS and 
HCO + imaging of these systems is valuable in confirming such a picture, a possible result of 
the different merger state of these two strongly interacting systems. 



3.3. Literature data 

A major goal of this study was to use our observations as well as results from the 
literature to compile a large, identical molecular line catalog for Arp 220 and NGC 6240. 
Thus the literature was scanned for all relevant molecular line measurements and for each 
transition the weighted average of every independent flux measurement was calculated. In 
cases where a measurement deviated significantly from the mean for obvious reasons (limited 
bandwidth, poor conditions, etc.), it was discarded from the average. Note that often in the 
literature line measurements are given only in (main beam or antenna) temperature units and 
only rarely in Janskys. In those cases, we used the aperture and/or main beam efficiencies 
quoted in the paper or, if those were not given, the average values for the particular telescope 
(typically tabulated in the observing manual). The final tally of lines is given in Table [3] 
along with their corresponding fluxes and global line luminosities. Global line ratios were 
calculated as the velocity/area averaged brightness temperature ratios (equivalent to taking 
the ratio of the line luminosities) and are listed in Tables H] and 

Our observations of 12 CO 2-1 in both Arp 220 and NGC 6240 yield somewhat higher line 
fluxes than previous single-dish and interferometry observations listed in Table [3], although, 
we note that Wiedner et al. (2002) using the JCMT measured an even larger line flux in 
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Arp 220 than ours. A careful scrutiny of the individual sub-scans in our data and of the 
calibration on the dates the observations were made revealed nothing out of the ordinary. 
In fact, given that each source was observed on different nights, yet both have ~ 40 per 
cent higher line fluxes than previously published values, suggest that the offset may well 
be real. A close inspection of previously published 12 CO 2-1 spectra reveals they all suffer 
from either limited bandwidth (the case for single-dish measurements), and/or the sources 
have been over-resolved spatially (in the case of interferometric observations). Because only 
few observations favor the higher 12 CO 2-1 fluxes in Arp 220 and NGC6240, however, we 
conservatively adopt the mean flux value of all available measurements for our analysis. The 
HCN 4-3 line in Arp 220 was first observed by Wiedner et al. (2002), who found a line flux 
only half as large as ours. Here, we adopt our own measurement since our spectrum has a 
higher signal-to-noise ratio than that of Wiedner et al. (2002), and the latter used aperture 
efficiencies which deviate significantly from the standard values given for JCMT. 

In NGC 6240 we find a HCN 1-0 flux that is half the previously published single-dish 
measurement (Solomon, Downes & Radford 1992a) but is in excellent agreement with the 
infererometric measurement by Nakanishi et al. (2005). Because our data were obtained over 
the course of two days in excellent weather conditions and because a careful reduction of the 
data from each day separately yielded consistent results, we find the lower values more likely. 
In the case of Arp 220, our estimate of the HCN 1-0 line flux is about 35 per cent larger 
than that of Solomon, Downes & Radford (1992a), which is just within the errors. Recent, 
independent observations of the HCN 1-0 line in Arp 220 and NGC 6240 with the IRAM 
30-m Telescope yield line fluxes in agreement with our measurements (Gracia-Carpio private 
communications), suggesting the HCN 1-0 line fluxes given by Solomon, Downes & Radford 
(1992a) may have been affected by systematic calibration uncertainties. Nonetheless, for the 
line ratio analysis in this paper, we conservatively adopt the weighted average of all three 
HCN 1-0 measurements in these two systems. 

We end this section by noting that Table [3] constitutes the most extensive molecular line 
catalog available for any (U)LIRG. In terms of extragalactic objects in general, it is matched 
only by data published for the much less extreme starburst galaxies M82 and NGC 253 (e.g., 
WeiB et al. 2005; Seaquist et al. 2006; Bayet et al. 2004; Bradford et al. 2003). However, 
unlike M82 and NGC 253, where only the 12 CO 1-0 and 2-1 transitions have been observed 
over their entire extent (e.g. WeiB, Walter & Scoville 2005), our study picks up the global 
emission from Arp 220 and NGC 6240 for all transitions, allowing for an unbiased study of 
the properties of their entire molecular gas reservoir using the corresponding mm/sub-mm 
line ratios. This mirrors observations of the mm/sub-mm lines observed in high- z starbursts 
where much fewer lines are available, and thus cannot yield meaningful constraints on the 
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conditions or the mass per gas phase for their molecular gas reservoirs. 

4. Analysis 

The first indication that the molecular ISM in extreme starburst galaxies is fundamen- 
tally different from that of normal galaxies came from the observation that the former have 
much larger 12 CO/ 13 CO line ratios (e.g. Casoli, Dupraz & Combes 1992). These are almost 
exclusively observed in ULIRGs and were interpreted as either due to missing 13 CO caused 
by a top-heavy initial mass function in these extreme systems (Casoli, Dupraz & Combes 
1992) or due to purely optical depth effects (i.e., low 12 CO optical depths) caused by a warm 
and highly turbulent H 2 medium present in these objects (Aalto et al. 1995). Then the bulk 
of the 12 CO emission arises from a diffuse enveloping phase where 12 CO has low to moderate 
optical depths (r ~ 1), and the large 12 CO/ 13 CO ratios are due to 13 CO-bright gas with a 
much a smaller filling factor than the diffuse 12 CO-bright phase. 

It is instructive to examine how far a single gas phase can go in terms of explaining the 
observed line ratios in Arp 220 and NGC 6240 before the fit becomes inadequate and a second 
gas phase must be introduced. To this end we employ a Large Velocity Gradient (LVG) model 
based on the work by Richardson (1985), which uses the observed CO line ratios to constrain 
the gas density (n(H 2 )), kinetic temperature (T fc ), and A co = [CO /H 2 }(dv / dr)' 1 , which is 
the CO abundance relative to H 2 divided by the cloud velocity gradient. The models span 
n(H 2 ) = 10 2 - 10 9 cm~ 3 (in steps of logn(H 2 ) = 0.5), T k = 10- 100 K (in steps of 5K), and 
A co = 10 6 — 10^ 3 (kms^pc" 1 )- 1 (in steps of logA co = 0.5). The model adopts the CO-H 2 
collisional rate constants from Flower & Launay (1985), and included all rotational levels up 
to J = 11 for kinetic temperatures of 10, 20, 40, 60, and 100 K were included. Collision rates 
for intermediate kinetic temperatures were estimated by interpolation. Finally, because the 
assumed 12 CO/ 13 CO abundance ratio can be important for whether a fit converges or not, 
we considered models with abundance ratios of 40, 60, and 100, values that span those found 
in the Milky Way (Langer & Penzias 1993). 

The need for a two-phase molecular ISM is perhaps most obvious in Arp 220, where 
no single LVG solution can account for even the 12 CO 1-0, 2-1, and 3-2 lines. In this 
system we have the interesting situation where the ratio r 21 = 12 CO (2-l)/(l-0)= 0.7 ± 0.1 
points towards sub-thermal excitation of even 12 CO 2-1 (n crit = 2.7 x 10 3 cm~ 3 ) where 
typically n(H 2 ) ~ 3 x 10 2 cm -3 and T k ~ 40 - 60 K fits such a ratio well (for LTE such a 
low ratio implies < 7K the minimum gas temperature permitted by cosmic-ray heating, 
Goldsmith & Langer 1978). Such conditions cannot reproduce the observed r%2 = 12 CO(3- 
2)/(l-0)= 1.0 ± 0.1 which indicates much more excited gas with density > 3 x 10 3 cm -3 and 
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kinetic temperature > 30 K. Going ahead to searcg for a LVG solution constrained by the 
12 CO 1-0, 2-1, 3-2 and just one 13 CO line strength, the J=l-0 lines, yields predictably a 
much worse solution given by n(H 2 ) — 1 — 3 x 10 2 cm -3 and ~ 35 — 45 K. Aside from 
poorly reproducing the aforementioned line ratios, barely satisfying the Tj. > T^ust criterion 
for these galaxies (see 4.2), such a gas phase corresponds to brightness temperature 13 CO 
(2-l)/(l-0), (3-2)/(l-0) ratios of r 2 i( 13 CO) ~ 1.4 - 1.5 and r 32 ( 13 CO) ~ 0.6 - 0.8, in sharp 
contrast to the observed ones of ~ 1.64 and ~ 5 respectively. 

In the case of NGC6240, both r 2 i = 1.2 ± 0.2 and r 32 = 1.1 ± 0.2 ratios indicate higly 
excited gas with a typical gas phase model of n(H 2 ) — 1 x 10 4 cm 3 and Tk — 40 — 50 K, that 
can also fit the observed large J=l-0 12 CO/ 13 CO ratio (7£io(model) = 37). Nevertheless 
this solution fails dramatically for the J=2-l, J=3-2 12 CO/ 13 CO line ratios where it yields 
IZ21 ~ ~ 17, much lower than the observed values (Table HJ), even for the largest 
plausible abundance of [ 12 CO/ 13 CO]= 100. This could not have been found without our 
sensitive 13 CO 2-1 and 3-2 measurements, which in the case of NGC6240, have pushed the 
7^21 and IZ32 ratios to the highest such values reported in the literature. 

Finally, using different [ 12 CO/ 13 CO] abundance ratios namely [ 12 CO/ 13 CO] = 40, 60, 100, 
encompasing the range deduced for the Milky Way (Langer & Penzias 1993; Casassus, Stahl 
Sz Wilson 2005) we still find no single set of average molecular gas physical conditions that 
comes even close in reproducing the 12 CO and the 13 CO line ratios in either galaxy. 

4.1. C ls O emission and the opacity of the 13 CO lines 

The emergence of 13 CO transitions as the first ones marking a strong deviation from 
any single gas phase reproducing the relative strengths of the 12 CO transitions is further 
highlighted in the case of Arp220 by the observed strong C 18 J=l-0 line emission. There 
the C 18 0/ 13 CO 1-0 line ratio of 1.0 ± 0.3 is the highest such ratio obtained over an entire 
galaxy — by comparison, Papadopoulos, Seaquist & Scoville (1996) measure such high values 
towards only a few of the brightest GMCs in NGC 1068. Assuming the low abundance of 
[C 18 0/ 13 CO] = 0.15 found for the Milky Way (Langer & Penzias 1993) such high C 18 0/ 13 CO 
brightness ratios argue for much of the 13 CO J=l-0 emission emerging from a phase where 
r io( 13 CO) > 1, sharply deviating from the optical thin transitions deduced from the LVG 
of only the 12 CO, 13 CO transitions described in the previous section. A much denser gas 
phase can naturally yield the significantly higher optical depths needed to explain the high 
C 18 0/ 13 CO J=l-0 ratio. 
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4.2. The properties of the dense molecular gas 

In this section we use the many transitions of HCN, CS, and HCO + in our catalog - 
probing a wide range of critical densities and excitation temperatures in a redundant manner 
(see Table - to constrain the properties of the dense gas phase. 

In order to model the radiative transfer of these lines, we constructed a modified version 
of the original LVG code. The modifications were straightforward since HCN, CS, and HCO + 
are simple linear molecules like CO. For HCN, CS, and HCO + , the collision rate constants 
provided by Green & Thaddeus (1974), Turner et al. (1992), and Flower (1999) were adopted. 
We explored a similar density range as we did for the CO lines, but limited the temperature 
range to Td < Tk < 120 K. The lower limit set by the dust temperature (Td), reflects the 
fact that the dust and gas are thermally de-coupled, with the FUV-induced photoelectric 
and turbulent gas heating (and its cooling via spectral line rather than continuum emission) 
setting Tk £ Td- Only in the densest, most FUV-shielded and quiescent regions of the 
GMCs can we expect Tk ~ T^. The upper limit is the typical equilibrium temperature of 
purely atomic, C + -cooled HI in the cold neutral medium (e.g. Wolfire et al. 2003). Based on 
detailed fits to their FIR/sub-mm SEDs (Dopita et al. 2005), we adopted a dust temperature 
of ~ 40 K for both Arp220 and NGC6240. Since the molecules have different abundances, 
the range in A is also different for a common velocity gradient (dv/dr). In our own Galaxy, 
typical abundances are [HCN/H 2 ] ~ 2 x 10 -8 (Bergin, Snell & Goldsmith 1996; Lahuis 
& van Dishoeck 2000), [CS/H 2 ] ~ 10~ 9 (Paglione et al. 1995; Shirley et al. 2003) and 
[HCO + /H2] = 8x 10~ 9 (Jansen 1995). In order to accommodate the substantial uncertainties 
associated with the abundances, the range of A used was 10 -10 — 10~ 7 , 10~ n — 10~ 8 , and 
10~ n - lO" 8 (kms^pc -1 ) for HCN, CS, and HCO + , respectively. These correspond to 
velocity-gradients in the range ~ 0.2 — 800 km s -1 pc _1 for the abundances quoted above. 

In the case of Arp 220, the observed HCN ratios are consistent with LVG solutions in the 
range n(H 2 ) — (0.3 — 1) X 10 6 cm~ 3 , and A HCN ~ (0.3 — 1) x 10~ 9 (kms _1 pc _1 ) , for which 
Thcn(i-o) — 3 — 7. From the fit of the CS lines we find n(H 2 ) — (0.3 — 3) x 10 6 cm -3 , and A cs ~ 
(0.01 - 1) x 10~ 8 (kms^pc -1 ) . Finally, the subthermal HCO + ratios allow for solutions 
in the range n(H 2 ) ~ (0.3 — 1) x 10 4 cm~ 3 , and A HCO + ~ (0.3 — 1) x 10~ 8 (kms _1 pc _1 ) . 
All the aforementioned solutions are for Tk > 40 K and their details are shown in Table 7. 
There it can be seen that while HCN, CS, and HCO + all trace what in broad terms can be 
characterized as dense gas (i.e., n(H 2 ) > 10 4 cm~ 3 ), with the density range of the CS and 
HCN solutions being very similar. Moreover we find the CS and HCN emission tracing gas 
with densities ~ 100 x higher than that traced by the HCO + lines. This simply reflects the 
high HCN, CS and the low HCO + line ratios measured in Arp 220, further reinforced by the 
~ 5 times lower critical densities of the HCO + transitions than e.g those of HCN at the same 
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rotational level (see Table 6). These results are certainly consistent with Galactic surveys 
of CS lines, which have shown they are unmistakable markers of very dense, high-mass star 
forming cores (Plume et al. 1997; Shirley et al. 2003). 

For NGC6240, the HCN LVG fits converge mostly over n(H 2 ) ~ (1 - 3) x 10 5 cm~ 3 , 
T k > 40 K, and A HCN ~ 3 x 10" 10 - 1 x 10" 8 (kms^pc" 1 ) ~\ with an optically thick HCN 
J=l-0 line (t hcn(1 _ 0) ~ 3 — 29). From the HCO + lines we find solutions in the range n(H 2 ) ~ 
(0.3 - 1) x 10 4 cm" 3 , T k > 40 K, and A HCO + - (0.3 - 1) x 10~ 8 (kms^pc -1 )" 1 - Thus in 
NGC6240, as in Arp220, HCN is tracing gas which is significantly denser than that traced 
by HCO + , over a similar range of temperatures. 

In an attempt to "break" the significant degeneracy that still exists in the LVG solutions 
for the dense gas we make use of the fact that at least in the Galaxy, most of the dense star- 
forming gas is found in self-gravitating cores (Shirley et al. 2003). In the cases of Arp220 
and NGC6240, we use this to try and limit the allowed (n(H 2 ), T k , A x )-parameter space by 
examining the ratio of the velocity gradient infered from the LVG fits to that expected for 
virialized cores, namely 

K _ (dv/dr) LVG [X/H 2 ] / n(H 2 ) V 1/2 

mr ~ {dv/dr) mr WEA X Vl0 3 cm^ ' W 

where [X/H 2 ] is the abundance ratio of the given molecule X (see Papadopoulos & Seaquist 
1999, or Goldsmith 2001 for a derivation), and a ~ 1 — 2.5 is a constant whose value depends 
on the assumed density profile of a typical cloud (Bryant & Scoville 1996). 

For virialized gas we expect K vir ~ 1 to within a factor of a few (due to uncertainties in 
cloud geometry, density profile, and abundance assumed to deduce the observed dv/dr from 
A x ). LVG solutions that correspond to K V i T » 1 indicate non-virial, unbound motions and 
are ruled out for dense star-forming gas along with those that have K V { r <C 1 (gas motions 
cannot be slower than those dictated by its self-gravity). For Arp220 and Table [7] we see 
that for both HCO + and HCN, LVG solutions with T k > 40 K have K vir ~ 1, and thus 
this conditions does not single out any particular range of solutions. For CS however, we 
find K V i r ^ 1 for most solutions except for T k = 35 — 40 K, n(H 2 ) = 3 x 10 6 cm -3 and 
A cs = 10~ 10 (kms -1 pc -1 ) -1 . This suggests that the true CS abundance may be significantly 
higher than assumed (by at least an order of magnitude) . Finally for NGC6240 all solutions 
that fit the HCO + ratios are consistent with K vir ~ 1, while for HCN this is the case for 
T k = 40 - 55 K and n(H 2 ) = 10 5 cm" 3 and A HCN = 10~ 8 (kms^pc^ 1 )" 1 (see Table 7). 
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4.2.1. HNC vs. HCN 

Observations of HCN and HNC in Galactic GMCs (Goldsmith et al. 1981), as well as 
steady state chemical models (Schilke et al. 1992), predict an increase in the HCN/HNC 
abundance ratio with increasing temperature and density, and attribute this to the destruc- 
tion of HNC (but not HCN) via neutral-neutral reactions with atomic oxygen or hydrogen. 
Then, at densities ~ 10 4-5 cm -1 , temperature-independent ion- neutral reactions, which form 
HCN and HNC in equal amounts, may become dominant (Schilke et al. 1992; Aalto et al. 
2002), and as a result the HNC abundance can be substantial. However, given the high 
optical depths implied for the HCN J=l-0 transition from our excitation analysis of the 
HCN line ratios, no useful limits can be placed on the [HNC/HCN] abundance ratio from 
the observed brightness temperarure ratio of 1.2 ± 0.2 measured for their J=l-0 transition 
in Arp220 (consistent with both lines being optically thick). For a comparative analysis of 
HCN versus HNC excitation using several rotational transitions, and the diagnostic nature 
of the [HNC/HCN] abundance ratio for the excitation environment of ULIRGs the interested 
reader is referred to recent work by Aalto et al. (2007). 

4.3. Sub-thermal HCO + line excitation: HCO + as dense gas mass tracer 

In both galaxies the low HCO + line ratios result in LVG solutions with densities mostly 
~ 10 — 100 times lower than those deduced from the analysis of the HCN ratios (Table [7]). 
Thus if only HCO + line observations were used, they would imply a bulk dense gas phase 
dominated by less extreme densities than is actually the case, especially in Arp220. The 
several HCO + and HCN transitions observed, and the independent analysis afforded by the 
CS transitions (for Arp220), make it unlikely that the aforementioned conclusions can be 
attributed to enhanced HCN abundances (and thus thermalization of its transitions at lower 
densities because of radiative trapping) caused by XDRs, or simply to IR-pumping of HCN 
J=l-0. Thus recent claims that HCO + is a better tracer of the dense gas in (U)LIRGs than 
HCN (Gracia-Carpio et al. 2006), do not seem to hold in the cases of Arp 220 and NGC 6240, 
and a more comprehensive approach is needed to decide the case (Papadopoulos 2007b). 

4.4. The mass of the dense gas phase in Arp 220 and NGC 6240 

In this section we make use of the constraints on the dense gas in Arp 220 and NGC 6240 
obtained in §4.21 to derive the total mass of the dense gas in these two systems. We do this 
using two different methods, the results of which are summarized in Table El The first 
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one assumes the dense gas reservoir to be fully reducible to an ensemble of self-gravitating 
units whose molecular line emission does not suffer any serious cloud-cloud shielding. Then 
it can be shown (e.g. Dickman, Snell, & Schloerb 1986) that the total mass of e.g. an 
HCN-luminous dense gas phase can be found from 



where a H cN — 

2.1 v /^)/T b , 

hcn(i-o) is a conversion factor (e.g. Radford, Solomon & Downes 
1991a). We insert the range of n(H 2 ) and TJ, iHCN(1 . 0) values inferred from the LVG solutions 
to the HCN line ratios, which for Arp 220 are n(H 2 ) = 0.3 x 10 6 cm -3 and T^HCNfi-o) = 
34 - 62 K, and for NGC6240 n(H 2 ) = (0.1 - 0.3) x 10 6 cm" 3 and T 6 , H cn(i-o) = 31 - 40 K. 
From these values we find: a HCN — (19 — 34) M® (Kkms -1 pc 2 ) -1 and a HC N — (17 — 
37) M Q (Kkms-ipc 2 )- 1 for Arp 220 and NGC6240, which in turn yields dense eas mass 
traced by HCN of M donsc (H 2 ) ~ (1.8 - 4.2) x 10 10 M & and ~ (1.0 - 2.8) x 1O 1O M for these 
two systems respectively. 

In a similar fashion, the analysis of HCO + and CS line ratios with their corresponding Tj, 
and n(H 2 ) values yields corresponding conversion factors which, along with the luminosities 
of the lowest HCO + and CS transitions available, yield the dense gas mass traced by each 
species (Table |8]). The constraints of an underlying density-size cloud hierarchy (see section 
4.5) demand M(H 2 ) HCN/CS < M HCO + (since the HCN, CS lines are tracing a denser gas than 
those of HCO+). Thus from Table 8 we obtain M dense (E 2 ) ~ 1-5 x 1O 1O M (Arp 220) and 
Mfen.se(H 2 ) ~ (1 — 2) x 1O 1O M (NGC6240) as the best dense gas mass estimates conforming 
to the aforementioned inequality. Systematic biases can be introduced by a) the lack of a 
pressure-term in eq. [6] correcting for molecular gas overlying the dense gas regions, b) by the 
existence of substantial stellar mass even within the dense gas phase. In both of these cases 
the true conversion factor and the deduced gas masses will be lower (see Bryant & Scoville 
1996; Downes & Solomon 1998 for the appropriate formalism). 

Finally, given the significant number of HCN transition available for both systems, we 
can estimate the dense gas mass by assuming the bulk of the HCN molecules are in rotational 
states J < 6. Then 



where Nj is the total number of HCN molecules in state J and jj = 1.36 accounts for He. 
The HCN population numbers are calculated from 
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M donsc (H 2 ) = [H 2 /HCN] [JVj + N 2 + N 3 + N, + N 5 + N 6 ) /xm : 
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where Aj 5 j_i and = (1 — e~ Tj -'- 1 ) / Tj t j-i are the Einstein coefficient and the escape 

probability for the J — > J — 1 transition, respectively. Unlike the first method this approach 
depends directly one the assumed [HCN/H 2 ] abundance ratio while it is in principle indepen- 
dent of any assumptions regarding the velocity field of the dense gas phase (i.e. it remains 
valid for a non-virialized, unbound gas phase), as long as the HCN line emission remains 
effectively optically thin (i.e. any large r's emerge within small gas "cells"). The dense gas 
masses estimated with the second method (see Table 8) bracket the best estimates yielded 
by the first one (constrained by the density- mass hierarchy assumtion), though more tightly 
for Arp 220 than NGC 6240. 

The two methods overlap over the same range of M(H 2 )dense ~ (1 — 2) x 1O 1O M for 
both galaxies, which is comparable to the total molecular gas mass estimated independently 
from models of their interferometrically imaged CO emission (Downes & Solomon 1998), or 
their CI J=l fluxes (Papadopoulos & Greve 2004). Thus the bulk of the molecular ISM in 
these two extreme starbursts is in a dense state with n(H 2 ) ~ (10 5 — 10 6 ) cm~ 3 . 

It is worth comparing this mass to the fundamental upper limit set by the dynamical 
mass in these objects. For Arp 220, the latter was estimated to be M dyn ~ 4 x 10 10 M Q 
(adopted to the cosmology used here) within a radius of 1.4 kpc, based on resolved obser- 
vations of CO (Downes & Solomon 1998). For NGC 6240 M dyn ~ (0.7 - 1.6) x 1O 1O M 
within a radius of ~ 500 pc (Tacconi et al. 1999; Bryant & Scoville 1999). The CO emit- 
ting region in NGC 6240 extends significantly beyond that (Tacconi et al. 1999) and thus a 
significant amount of molecular gas may reside beyond the inner 500 pc. Assuming a con- 
stant rotation curve, we obtain an upper limit of M dyn = 6 x 10 10 M Q within R = 1.4 kpc 
for NGC 6240. Since any HCN/CS/HCO+ bright region will be contained within the CO- 
bright regions we expect M dyn > M co > M dense in both systems. Thus the range of 
M(H. 2 ) de nse ~ (1 - 2) x 10 w M Q corresponds to ~ 20% - 50% of the dynamical mass in 
Arp 220, and > 15% up to 100% for NGC 6240. 



4.5. The density-size relation for the dense gas: steeper than Larson's law? 

In this section we attempt to go a step beyond using HCO + , HCN, and CS lines as 
isolated estimators of the dense gas mass by assuming their emission emerging from an 
underlying density-size hierarchy, found for Galactic molecular clouds by Larson (1981), and 
verified by subsequent high-resolution multi-transition studies. The origin of the density-size 
and linewidth-size power laws revealed for GMCs in the Galaxy may lie in the characteristics 
of supersonic turbulence and its global driving mechanisms (Heyer & Brunt 2004) when 
present on self-gravitating structures. As such they may remain similar even in the extreme 
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ISM conditions present in ULIRGs. 

A rigorous approach to test this in extragalactic environments (where molecular clouds 
cannot be viewed with the same level of detail) using the global molecular line emission of 
various species would require incorporating such power laws in the radiative transfer model 
used to interpret the emergent molecular line emission. In our case we opt for a simpler 
approach in which the denser HCN-bright gas is simply assumed to be "nested" within the 
less dense regions traced by the HCO + transitions. 

We can rudimentary test for this by deriving the velocity-averaged area filling factor of 
the HCN-bright relative to the HCO + -bright gas phase: / H cn,hco+, which ought to be < 1. 
This can be found from 

T ' T 

HCN _ f i fc,HCN / Q x 

t, — J HCN, HCO+ rp ■ \ y ) 

L HCO+ J6,HCO+ 

For the LVG solution ranges above T k = 40 K (Table 7) and the J=l-0 transition of both 
species, the observed line luminosities and LVG-deduced intrinsic brightness temperatures 
yield: 0.2 < / HCN3CO+ < 1.0 (Arp220) and 0.1 < / HC n,hco+ £ 0.3 (NGC6240). For an under- 
lying density-size power law (n(H 2 )) oc R~ a , we set a = 2log (n H co+/ n HCN) /logf ilCN]HCO + as 
an approximation (n HCN , n HCO + are the mean volume densities deduced from the LVG fits 
of the HCN and HCO + line ratios). For the range of densities corresponding to LVG fits 
of the HCN and HCO + ratios (for Tk > 40 K), we find a > 2, steeper than the value of 
a ~ 1 typical for Galactic molecular clouds. Using the HCO + and HCN J=4-3 transitions 
instead yields 0.09 < / H cn,hco+ £ 0.6 (Arp220), and 0.03 < / HCN3CO+ < 0.2 (NGC6240). These 
smaller values are expected for emission from the two lines bracketing the entire range criti- 
cal densities in our line inventory (Table 6), and an underlying density-size ISM hierarchical 
structure. However, even for these values of / H cn,hco+ the corresponding density ranges yield 
a > 1.3, still steeper than a Larson-type power law. Such a steepening could in principle 
be the result of extreme tidal stripping of GMCs and then merging of the surviving densest 
clumps to super-dense structures, yet another signature of vastly different average molecular 
gas properties in extreme starbursts than those in more quiescent environments. 



4.6. The diffuse gas phase 

The observed 12 CO line emission will have contributions from both the dense and diffuse 

gas and the observed global line ratio of any given J + 1 — > J line with respect to 12 CO 
1-0 can therefore be written as 

(A) r (12) (B) 

'j+l,J~r u 10 ' J+1,J / in x 

r w = : . ^(12) ' ( 10 ) 
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where rj^i j is the line ratio for the diffuse phase (denoted by A), rf^ x 3 the one for the dense 
phase (denoted by B), and the contrast factor, C\ is discussed below. It is seen from eq. [10] 
that in order to constrain the properties of the diffuse gas, one has to determine and subtract 
the contribution from the dense gas phase. Fortunately, this can now be done with the 
constraints put on the dense phase properties in §4.21 The contrast factor, which determines 
the relative emission contribution of the two gas phases, is given by = /ba^^/T^), 
where /ba is the velocity-integrated relative filling factor between the two phases and 
are the velocity- and area-integrated line brightness temperatures of 12 CO 1-0. In order 
to determine this quantity, we make the reasonable assumption that the HCN emission is 
dominated by the dense gas phase and, as a result, we can write the global (observed) 
HCN(1-0)/ 12 CO(1-0) line ratio as 

rp r {12) T (B) 

J fe,HCN _ ^lO b, HCN 

We then estimate using the observed value of T bjHCN /T bjCO and the ^hcn/^&co va l ues 
obtained from the LVG fits of the dense gas phase. The observed global line ratio for Arp 220 
is Tfe iHCN /Tfe iCO = 0.18 ± 0.03, and ^ hcn/^ co = 0-5 ~ 0-8 as estimated from the dense gas 
solution range. In the case of NGC 6240 the same two line ratios are 0.08 ±0.01 and 0.4 — 0.8. 
Folding in the uncertainties on the observed global ratios and allowing for the range in the 
ratio predicted for the dense gas, we find C\ Q ~ 1.5 — 3.0 (Arp 220) and C{ ~ 1.3 — 2.8 
(NGC 6240). 

With Cffl determined, we can then disentangle the contributions to the 12 CO line ratios 
from the diffuse and dense phases using Eqs [101 If we assume the TZ ratios dominated by the 
diffuse phase, the resulting 12 CO and 13 CO emission line ratios for Arp 220 are consistent 
with n(R 2 ) ~ (1 - 3) x 10 2 cm " 3 , T k > 40 K, and A co ~ (3 - 10) x 10" 5 (km s' 1 pc -1 ) -1 . 
For NGC 6240, the LVG model converges on n(H 2 ) ~ (1 - 3) x 10 3 cm~ 3 T k > 40 K, and 
A co ~ (1 — 3) x 10~ 6 (kms -1 pc -1 ) -1 . All solutions correspond to non-virialized K vir > 5 
values (assuming a Galactic abundance [CO/H2]= 10~ 4 ) and have moderate 12 CO J=l 
optical depths (r 10 ~ 1). This seems to be the typical diffuse gas phase found in such systems 
(e.g. Aalto et al. 1995; Downes & Solomon 1998) which, while dominating their low-J CO 
line emission, it does not contain much of their total molecular gas mass. 



5. Dense gas mass, star formation, and its efficiency in Arp 220 and NGC 6240 



The dense gas masses derived for Arp 220 and NGC 6240 in §4.41 allow us to derive accu- 
rate star formation efficiencies in these two galaxies. Although the star formation efficiency 
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(SFE) was originally defined as SFE = SFR/M(H 2 ), where SFR is the star formation rate 
and M(H 2 ) is the total molecular gas mass, we shall here focus on the SFE per dense gas 
mass, with the SFR parameterized by the IR luminosity (i.e. SFE oc L m /M dense ). 

Using the original definition of IR luminosity as the SED integrated from 40- to 120- 
Hm (Helou et al. 1986), we derive L m = 1.6 x 1O 12 L (Arp220) and L m = 7.0 x 10 11 L Q 
(NGC6240). For the deduced dense gas mass range of M dense (R 2 ) ~ (1 - 2) x 10 10 M we 
obtain L m /M dense ~ (80 - 160) L & M & ~ l (Arp220) and L m /M dense = (35 - 70) L Q M Q " x 
(NGC6240), suggesting a somewhat higher star formation efficiency in Arp220. Finally 
these values are within the range of those found for a large sample of infra-red galaxies, 
consisting of large spiral galaxies, LIRGs, and a few ULIRG by Gao & Solomon (2004a,b) 
using the HCN 1-0 line as a linear tracer of dense gas. 

A density-size hierarchy found for Galactic GMCs remaining valid in other galaxies 
would suggest that HCO + , HCN, and CS lines, given their increasing excitation require- 
ments, would probe progressively deeper into the warm and dense gas phase fueling star 
formation. We briefly explore this scenario by calculating three separate star formation 
efficiencies, SFE HCO +, SFE HCN , and SFE CS , using the separate dense gas mass estimates 
based on HCO + , HCN, and CS available for Arp220 only. For this galaxy we find SFE CS 
(~ 107-533 L Q M q _1 ) is larger than SFE HCN (~ 40-94 L & M _1 ), which may suggest that 
the star formation increases as we probe higher densities (Shirley et al. 2003). On the other 
hand, this trend is not confirmed by HCO + , which yields SFE HCO + ~ 100 — 320 L Q M Q -1 , 
just as we do not see a significantly narrower range in SFE-values at higher densities as 
one might otherwise expect (Shirley et al. 2003). In NGC6240 the star formation efficiency 
derived from HCO+ (SFE HCO + ~ 32 — 54 L & M & x ) is smaller than the efficiency obtained 
from HCN (SFE HCN ~ 41 - 70 L Q M Q _1 ). We conclude that given the uncertainties in the 
mass estimates we cannot robustly claim a significant increase in the SFEs derived using 
HCO+, HCN and CS. 

Theory currently provides some hints that support a small range of star formation 
efficiencies per dense gas mass. Simulations of star formation in turbulent media where 
dense (n(H 2 ) > 10 4 cm~ 3 ) gas clumps decouple from a supersonic, turbulent cascade to form 
stars at a constant efficiency seem to reproduce well some of the observed characteristics of 
SFR versus gas in galaxies (Krumholz & McKee 2005). Furthermore, Scoville (2003) showed 
that the negative feedback effect of high-mass star formation naturally leads to an upper 
bound on the SFE for star forming regions of 500 L & Mq -1 . Although the exact value 
of this upper limit may change in the future (since it depends on the poorly constrained 
effective radiative opacity of the dust), the underlying physics is expected to apply to local 
star forming regions as well as those in circumnuclear starbursts. We note that our derived 
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SFEs are consistent with the current upper limit. 

6. Implications for dense star forming gas at high redshifts 

Evidence from Galactic observations and recent results from the ULIRG/QSO Mrk231 
(Papadopoulos, Isaak & van der Werf 2007) point to the same gas phase responsible for both 
HCN and high-J CO line emission. This is expected, as such a phase is intimately linked 
with ongoing star formation and thus its warm and dense conditions will excite the high-J 
CO transition and those of HCN, CS and HCO + alike. Thus the weaker lines of heavy 
rotor molecules with frequencies within the numerous transparent atmospheric windows at 
< 350 GHz, can be used to probe the gas emitting also the luminous CO J+l— > J, J+ 1 > 4 
lines at frequencies of v > 460 GHz beyond which the atmospheric transparency window 
rapidly closes, even in sites where prime sub-mm telescopes are located. 

In Fig. [3](a) and (c) we have plotted the expected line fluxes of the 12 CO and 13 CO line 
templates - based on Arp220 and NGC6240, respectively - for redshifts z = 0.05,0.1,0.5, 
and 1.0. The 12 CO and 13 CO line fluxes were extracted from the LVG model with input 
conditions set by the HCN solution ranges for the dense gas in these two systems (Table [7]). 
Similarly, Fig. EIb) and (d) show the HCN/CS/HCO + line templates for the solution ranges 
given in Table [7J For comparison we also show the latest sensitivity estimates for the HI-FI 
instrument (de Graauw et al. 1998, 2005) - a series of highly advanced heterodyne receivers 
due to fly onboard Herschel, an ESA cornerstone mission expected to launch in 2008. HI-FI 
will observe in six bands covering the frequency range 480-1910 GHz, enabling it to target 
the highest J-transitions in \ow-z sources. Fig.[3]shows that if Arp 220 and NGC 6240 are not 
significantly different from the bulk population of (U)LIRGs in the redshift range < z < 1, 
HI-FI/Herschel will be able to detect the CO J — > J — 1, J > 5 lines in such objects out to 
z ~ 0.1 in less than 1 hr. At z ~ 0.5 the CO J —>■ J — 1, J > 7 lines are still within the HI-FI 
frequency range, and are easily detectable. We see the HCN and HCO + J — ► J — 1, J > 6 
lines are detectable in systems akin to Arp 220 and NGC 6240 out to z ~ 0.1. At z ~ 0.5, 
however, HCN lines are still detectable only in Arp220-like systems. The predicted HCO + 
lines fluxes for both type of systems at z ~ 0.5 fall below the 1 hr 5-cr sensitivity limits of 
HI-FI. 

Thus for (U)LIRGs in the redshift range < z < 0.5, HI-FI will complement ground- 
based (sub)mm facilities where only the low- J lines are accessible and thereby densely sample 
the rotational ladder for a number of important molecules in such systems. This will allow 
a full inventory of the molecular ISM in luminous starbursts and (U)LIRGs and provide a 
valuable set of low-z templates for interpretation of high-J detections in extremely distant 
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galaxies (z > 2) with ALMA. Characterizing these transitions relative to the lower ones in 
the local extreme starburst population - similar to the need to understand local dust SEDs 
before interpreting the very few frequency points detected for the dust continuum in high- 
z galaxies - will prove to be extremely valuable before embarking on the high- J, high- ,2 
ventures with ALMA. 



7. Summary 

The main results presented in this paper are: 

1. Using the JCMT and IRAM 30-m telescopes we report first detections of a number of 
important molecular lines, including transitions of 13 CO, HCN, CS, and HCO + , in the two 
prototypical ULIRGs Arp220 and NGC6240. These observations along with measurements 
from the literature complete what is currently the largest molecular line catalog for such 
systems, capable of probing a wide range of densities and temperatures, especially for their 
important dense and warm star-forming gas phase. 

2. In Arp220 a systematic change in line profile as higher densities are probed suggests 
large scale gas phase segregations within this ULIRG. This is not observed in NGC6240 
even in the profiles of lines widely separated in their excitation requirements such as 12 CO 
J— 2-1 and HCN J= 4-3. This points to the exciting possibility of using high-J/high-dipole 
molecular lines to deduce the state of gas relaxation in ULIRGs by simply comparing the 
profiles of lines with widely different excitation requirements. This could prove important 
since currently only in the closest ULIRGs can interferometry be relied upon to reveal the 
state of the gas directly. 

3. A single molecular gas phase is inadequate to account for the observed properties 
of the molecular ISM in Arp220 and NGC6240. Analysis of the numerous line ratios using 
a two-phase LVG model finds a diffuse (n(H 2 ) ~ 10 2_3 cm~ 3 ) phase that is most likely un- 
bound, and a dense gas phase (n(H 2 ) ~ 10 5 ~ 6 cm~ 3 ) that contains the bulk of the molecular 
gas in these systems {M(B. 2 ) dense ~ (1 - 2) x 10 10 M ). 

4. The four highly excited HCN rotational transitions (also four CS transitions for 
Arp220) correspond to a gas phase that is ~ 10 — 100 times denser than that suggested 
by the low, sub-thermal, HCO + line ratios measured in both systems. Thus HCN trace a 
denser gas phase than the HCO + lines. Moreover, these two phases are found consistent 



-24- 



with an underlying density-size power law, albeit with a steeper exponent than that found 
in Galactic GMCs (i.e. the Larson's law for densities). We argue that this could be another 
result of the extreme tidal stripping of GMCs in these strongly interacting systems. We 
also find that for CS all realistic solutions imply K vir << 1, which suggest that the true CS 
abundance is significantly higher (by an order of magnitude) than the assumed abundance 

(\cs/h 2 ] = 10- 9 ;. 

5. Finally, using Arp220 and NGC6240 as templates for starburst galaxies at higher 
redshifts, we show that the HI-FI instrument onboard Herschel will be ideal for detect very 
high-J transitions of CO/HCN/CS/HCO + in such objects out to redshifts of z ~ 0.5 while 
still preserving the all important information about their relative velocity profiles. Thus the 
Herschel Space Observatory, in conjuction with an ongoing large mm/sub-mm molecular line 
survey of extreme starbursts, will be able to produce the first densely sampled rotational 
ladders of CO, HCN, and HCO + for (U)LIRGs in the local Universe, a prerequisite in 
understanding star formation and its relation to dense gas in extreme starbursts localy and 
in their high redshift counterparts. 

We thank the anonymous referee for a thorough report which helped improve the paper 
substantially. We are grateful to all the telescope operators at the IRAM 30-m and the long- 
standing support we enjoy from the JCMT support personel and telescope operators. We are 
particular thankful to Jim Hoge at the JCMT for his active role in getting the observations 
done, often during adverse technical conditions. 
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Fig. 1. — Molecular line spectra observed in Arp220. The 12 CO 2-1 and 3-2 spectra have 
been scaled down by a factor six in order to fit on the plot. The velocity- axis is relative to 
^lsr — 5454 km s -1 . The FWHM line width and channel-to-channel rms noise are given for 
each line. All spectra were obtained with the JCMT except for those in panel c), e), h)-i), 
and m)-n) which are from the IRAM 30-m Telescope. 
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Fig. 2. — Molecular line spectra observed in NGC6240. The velocity-axis is relative to 
^lsr = 7359 km s -1 , and the 12 CO 2-1 and 3-2 spectra have been scaled down by a factor 
eight. All spectra were obtained with the JCMT except for those in panel c) and e) which 
are from the IRAM 30-m Telescope. Other details are the same as in Fig. [H 
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Fig. 3. — Top panels: integrated line flux strengths of the full 12 CO/ 13 CO (left) and 
HCN/CS/HCO + (right) rotational ladders predicted by the dense gas template of Arp220 
as given by the solution ranges in Table [71 The line fluxes are shown for Arp 220-like systems 
at redshifts z = 0.05,0.1, and 0.5. In the case of 12 CO and 13 CO we have also shown the 
z — 1.0 template. Bottom panels: Same as above, except the constraints on the dense gas 
in NGC 6240 have been used as a template. Horizontal bars are the expected 5-<r sensitivity 
limits of the five HI-FI bands after 1 hr of integration. In panels a) and c), the red and blue 
curves correspond to the line templates for 12 CO and 13 CO, respectively, while in panels 
b) and d), they correspond to HCN and HCO + , respectively. The green curve in panel b) 
correspond to the line template for CS. 
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Table 1. Observational parameters for the lines observed with the JCMT. 





12 CO(2-l) 


12 CO(3-2) 


13 CO(2-l) 


13 CO(3-2) 


HCN(3-2) 


HCN(4-3) 


HCO+(4-3) 


CS(5-4) 


CS(7-6) 


Arp220 




















v obs [GHz] 


226.439 


339.648 


216.481 


324.711 


261.159 


348.203 


350.392 


240.581 


336.787 


Tsys [K] 


~430 


350-500 


390-580 


1020-1980 


500-600 


250-360 


260-270 


280-340 


300-330 


tint H 


600 


600 


10800 


2400 


6000 


5400 


4800 


5400 


4200 


NGC6240 




















u oh3 [GHz] 


225.025 


337.527 


215.128 


322.682 


259.528 


346.027 


348.203 




334.683 


Tsys [K] 


~410 


~290 


400-480 


700-810 


280-290 


370-750 


310-460 




280-300 


tint M 


600 


600. 


13200 


3600 


3000. 


7800 


5400 




3000 


HPBW ["] 


20.8 


14.0 


21.3 


14.6 


18.4 


13.5 


13.5 


20.0 


14.0 


Vap 


0.60(0.58)° 


0.53(0.51)° 


0.60(0.58)° 


0.53(0.51)° 


0.54(0.54)° 


0.53(0.51)° 


0.53(0.51)° 


0.56(0.57)° 


0.53(0.51)° 



a) Standard JCMT aperture efficiencies relevant for the observing periods in question. Values in parentheses arc our own measurements. 



Table 2. Observational parameters for the lines observed with the IRAM 30-m. 





13 CO(1-0) 


13 CO(2-l) 


HCN(l-O) 


CS(2-1) 


CS(3-2) 


HNC(l-O) 


C 18 O(l-0) 


Arp 220 
















v obs [GHz] 


108.242 


216.481 


87.055 


96.238 


144.356 


89.051 


107.830 


Tsys [K] 


~430 


225-500 


95-125 


200-300 


200-300 


200-300 


200-300 


Unt H 


600 


4320 


4320 










NGC6240 
















^6, [GHz] 




215.128 


86.512 










Tsys [K] 




225-500 


96-125 










tint M 




7200 


7200 










HPBW ["] 


22.4 


11.4 


28.3 


25.1 


18.0 


27.1 


22.4 


Vap 


0.56 


0.41 


0.62 


0.59 


0.52 


0.59 


0.56 



Table 3. Molecular and atomic transitions detected towards Arp220 and NGC6240. 



Line 


FWHM" 


Arp220 

S v Av*>) L' c ) 


Telescope 


Ref. 


FWHM™ 


NGC 6240 
S^Av^ L' c ) 


Telescope 


Ref. 


12 CO(1-0) 


504 
480 


410 ±41 
384±115 d 
496 ± 99 
419 ± 36 


5.9 ±0.6 
5.5 ± 1.7 
7.2 ±1.4 
6.0 ±0.5 


OVRO 
OVRO 
IRAM 30m 


[1] 
[2] 
[3] 


370 
469 


314 ± 64 
324 ± 33 
322 ± 29 


8.3 ± 1.7 
8.6 ±0.9 
8.5 ±0.8 


IRAM 30m 
OVRO 


[3] 
[18] 


12 CO(2-l) 


447 

380 
420 


1071 ± 107 
1088 ± 330 d 
1118 ± 112 
1100 ± 220 
1730 ± 519 
1549 ±311 
1127 ±69 


3.9 ±0.4 
3.9 ± 1.2 
4.0 ±0.4 

4.0 ±0.8 
6.2 ± 1.9 
5.6 ± 1.1 

4.1 ±0.3 


OVRO 
IRAM 30m 
NRAO 12m 
OVRO 
JCMT 
JCMT 


[1] 
[4] 
[5] 
[6] 
[9] 

This work 


300 
410 


1220 ± 366 d 
1741 ± 350 
1492 ± 253 


8.1 ± 2.4 
11.6 ±2.3 
9.9 ± 1.7 


IRAM PdBI 
JCMT 


[19] 

This work 


12 CO (3-2) 


306 
510 


4094 ± 758 
4154 ± 975 
3700±1110 d 
3168 ± 634 
3674 ± 405 


6.6 ± 1.2 

6.7 ± 1.6 
5.9 ± 1.8 
5.1 ± 1.0 
5.9 ±0.6 


JCMT 
HHT 
JCMT 
JCMT 


[7] 
[8] 
[9] 

This work 


390 


3205 ± 642 


9.5 ± 1.9 


JCMT 


This work 


13 CO(1-0) 


369 


< 9.4 
9 ± 2 


< 0.15 


SEST 

IRAM 30m 


[5] 

This work 




6.5 ± 1.9 


0.19 ± 0.06 


IRAM 30m 


[20] 


13 CO(2-l) 


450 
390 


56 ± 17 d 
54 ± 11 
60 ± 13 
70 ± 16 
59 ±7 


0.22 ±0.07 
0.21 ±0.04 
0.24 ±0.05 
0.28 ±0.06 
0.23 ±0.03 


IRAM 30m 

IRAM 30m 
JCMT 


[10] 

ri 1 1 
[J-iJ 

This work 
This work 


360 


28 ±7 

25 ±5 

26 ±4 


0.20 ±0.05 
0.18 ± 0.04 
0.19 ± 0.03 


JCMT 

TT3 A 1\/T Qn.v, 

IrtAlvl oUlTl 


This work 
This work 


13 CO (3-2) 


410 


404 ± 86 


0.71 ±0.15 


JCMT 


This work 




< 79 


< 0.3 


JCMT 


This work 


12 C 18 O(l-0) 


280 


9 ± 2 


0.14 ±0.03 


IRAM 30m 


This work 












HNCO(5-4) 


280 


5 ± 2 


0.08 ±0.03 


IRAM 30m 


This work 












HCN(l-O) 


550 


35±ll d 
59 ± 12 


0.85 ±0.27 
1.44 ±0.29 


IRAM PdBI 
IRAM 30m 


[12] 

This work 




26±8 d 
14±4 d 


1.17 ±0.36 
0.63 ± 0.18 


IRAM 30m 
NMA 


[13] 
[21] 



Table 3 — Continued 



Line 


FWHM" 


Arp220 

S v Av b ~> L /c > 


Telescope 


Ref. 


FWHM" 


NGC 6240 
S„ Ad 6 ) L' c 1 


Telescope 


Ref. 




530 


38 ± 12 
48 ± 14 
44 ± 6 


0.93 ±0.29 
1.17 ±0.34 
1.08 ± 0.15 


IRAM 30m 
IRAM 30m 


[13] 
[22] 


340 
410 


13 ±3 
16 ±5 
15 ± 2 


0.58 ± 0.13 
0.72 ± 0.22 
0.66 ± 0.09 


IRAM 30m 
IRAM 30m 


This work 
[22] 


HCN(2-1) 


530 


144 ± 43 


0.88 ± 0.26 


IRAM 30m 


[22] 


410 


45 ± 13 


0.51 ± 0.15 


IRAM 30m 


[22] 


HCN(3-2) 


540 
530 


361 ± 73 
258 ± 77 
312 ± 53 


0.98 ±0.20 
0.70 ±0.21 
0.85 ± 0.14 


JCMT 
IRAM 30m 


This work 
[22] 


430 
410 


160 ± 33 
74 ±22 
160 ± 33 


0.80 ±0.16 
0.37 ±0.11 
0.80 ± 0.16 


JCMT 
IRAM 30m 


This work 
[22] 


HCN(4-3) 


590 


587 ± 118 


0.90 ± 0.18 


JCMT 


This work 


320 


138 ± 29 


0.39 ± 0.08 


JCMT 


This work 


HNC(l-O) 


516 ± 38 
407 


55 ± 11 
34 ±7 
43 ±8 
40 ±5 


1.28 ± 0.26 
0.79 ±0.16 
0.96 ±0.19 
0.94 ±0.11 


IRAM 30m 
SEST 
IRAM 30m 


[14] 
[15] 

This work 




14 ±7 


0.60 ± 0.30 


SEST 


[15] 


CS(2-1) 


339 


12 ±3 


0.24 ±0.06 


IRAM 30m 


This work 












CS(3-2) 


455 


22 ± 7 d 
26 ±5 
25 ±4 


0.20 ± 0.06 
0.24 ±0.04 
0.23 ±0.03 


IRAM 30m 
IRAM 30m 


[10] 

This work 












CS(5-4) 


460 


68 ± 15 


0.22 ±0.05 


JCMT 


This work 












CS(7-6) 


390 


111 ±24 


0.18 ±0.04 


JCMT 


This work 


240 


60 ± 14 


0.18 ± 0.04 


JCMT 


This work 


HCO+(1-0) 




20±6 d 

19 ±6 
22 ± 7 d > e 

20 ±4 


0.48 ±0.14 
0.46 ±0.14 
0.53 ±0.17 
0.49 ±0.09 


IRAM PdBI 
IRAM 30m 
IRAM 30m 


[12] 
[13] 
[16] 


420 


21 ±3 
25 ± 8 d ' e 
21 ±3 


0.93 ± 0.13 
1.11 ± 0.35 
0.95 ± 0.12 


NMA 

IRAM 30m 


[21] 
[16] 


HCO+(3-2) 




49 ± 15 d - e 


0.13 ±0.04 


IRAM 30m 


[16] 




47 ± 14 d - e 


0.23 ± 0.07 


IRAM 30m 


[16] 


HCO+(4-3) 


280 


106 ± 23 


0.16 ±0.03 


JCMT 


This work 


400 


71 ± 17 


0.20 ±0.05 


JCMT 


This work 



Table 3 — Continued 





Arp 220 






NGC 6240 






Line fwhm" 




Telescope 


Ref. FWHM° 


S v Av bS > L' c 1 


Telescope 


Ref. 


Cl( 3 Pi - 3 P„) 


1160 ±350 


JCMT 


[17] | 


600 ± 180 


JCMT 


[17] 



a) Line widths are in kms -1 . 

b) Line flux densities are in units of Jykms 1 . 

c) Line luminosities are in units of xlO 9 Kkms~ 1 pc 2 . 

d) We assume a 30 per cent error on the flux measurement. 

e) The fluxes are derived from the luminosities quoted by Gracia-Carpio et al. (2006). 

[1] Scovillc ct al. (1997); [2] Scoville et al. (1991); [3] Solomon ct al. (1997); [4] Radford ct al. (1991a); [5] Aalto ct al. 
(1995); [6] Dowries & Solomon (1998); [7] Liscnfeld et al. (1996); [8] Mauersberger et al. (1999); [9] Wiedner et al. (2002); 
[10] Solomon ct al. (1990); [11] Papadopoulos & Scaquist (1998); [12] Radford et al. (1991b); [13] Solomon ct al. (1992a); [14] 
Huttemeister et al. (1995); [15] Aalto ct al. (2002); [16] Graci'a-Carpio ct al. (2006); [17] Papadopoulos & Greve (2004); [18] 
Bryant et al. (1999); [19] Tacconi et al. (1999); [20] Casoli, Dupraz & Combes (1992); [21] Nakanishi et al. (2005); [22] Krips 
ct al. (2007). 
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Table 4. CO line ratios for Arp 220 and NGC 6240. 



Galaxy 


12 CO(2 - 1) 
12 CO(l - 0) 

»*21 


12 CO(3 - 2) 
i2CO(l - 0) 


12 CO(l - 0) 
i3CO(l - 0) 
K 10 


12 CO(2 - 1) 
13 CO(2 - 1) 


12 CO(3 - 
i3CO(3 - 

72-32 


2) 
2) 


C 18 O(l-0) 
i3CO(l - 0) 
77(18) 


Arp 220 
NGC 6240 


0.7±0.1 
1.2 ±0.2 


1.0 ±0.1 

1.1 ± 0.2 


43 ± 10 
45 ± 15 


18 ±3 
52 ±12 


8 ± 2 
> 32 




1.0 ±0.3 



Table 5. Observed and modeled HCN, CS and HCO+ line ratios for Arp 220 and 

NGC 6240. 



Galaxy 


HCN(2 - 1) 


HCN(3 - 2) 


HCN(4 - 3) 


HCN(1 - 0) 


CS(3 - 2) 


CS(5 - 4) 


CS(7 - 6) 


HCO+(3 - 2) 


HCO+(4 - 3) 


HCN(1 - 0) 


HCN(1 - 0) 


HCN(1 - 0) 


HNC(1 - 0) 


CS(2 - 1) 


CS(2 - 1) 


CS(2 - 1) 


HCO+(l - 0) 


HCO+(l - 0) 


Arp 220 




















observed 


0.8 ± 0.3 


0.8 ± 0.2 


0.8 ± 0.2 


1.2 ± 0.2 


1.0 ± 0.3 


1.0 ± 0.3 


0.8 ± 0.3 


0.27 ± 0.10 


0.33 ± 0.09 


modeled 


1.0 


0.9 


0.8 




1.0 


0.9 - 1.0 


0.8 - 0.9 


0.3 - 0.4 


0.4 


NGC 6240 




















observed 


0.8 ± 0.3 


1.2 ± 0.3 


0.6 ± 0.1 


1.1 ± 0.6 








0.24 ± 0.08 


0.21 ± 0.06 


modeled 


0.9 - 1.0 


0.8 - 0.9 


0.6 - 0.7 










0.3 - 0.4 


0.4 
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Table 6. Line excitation characteristics (from Jansen 1995). 



Transition 


''rest 


E u /k B 


n crit 




[GHz] 


[K] 


[cm" 3 ] 


12 CO(1-0) 


115.271 


5.5 


4.1 x 10 2 


12 CO(2-l) 


230.538 


16.6 


2.7 x 10 3 


12 CO(3-2) 


345.796 


33.2 


8.4 x 10 3 


HCN(l-O) 


88.632 


4.3 


2.3 X 10 5 


HCN(3-2) 


265.886 


25.5 


4.1 x 10 6 


HCN(4-3) 


354.734 


42.5 


8.5 x 10 s 


HCO+(1-0) 


89.189 


4.3 


3.4 x 10 4 


HCO+(3-2) 


267.558 


25.7 


7.8 x 10 5 


HCO+(4-3) 


356.734 


42.8 


1.8 x 10 6 


CS(2-1) 


97.980 


7.1 


8.0 X 10 4 


CS(3-2) 


146.969 


14.1 


2.5 X 10 5 


CS(5-4) 


244.936 


35.3 


1.1 x 10 6 


CS(7-6) 


342.883 


65.8 


2.9 x 10 6 



"crit = A u i I ^2,i Cui i where C u i is the colli- 
sional rate 'out' of the level u (cither up or down), 
calculated for Tj, — 100 K in the optically thin limit. 



Table 7. LVG solution ranges, see §4.21 for details. The most likely solutions have been 

highlighted in bold. 





T fc /[K] 


n(H 2 )/cm- 3 


A/(km s 1 pc 


I)-! 


xV 


K vir 




Arp220: 














HCO + 
















10 - 40 


(1 - 3) x 10 4 


(0.3 - 1) X 10~ 


8 


2.1 - 2.8 


0.4 - 0.7 




45 - 95 


1 X 10 4 


0.3 X 10" 8 




2.1 


1.2 




100 - 120 


0.3 X 10 4 


1 X 10~ 8 




2.0 - 2.1 


0.7 


HON 
















10 - 25 


(0.3 - 1) X 10 6 


(0.3 - 1) X 10- 


7 


0.5 - 1.3 


0.02 - 0.17 




30 - 40 


1 X 10 6 


0.3 X 10~ 9 




0.7 - 0.8 


3.1 




45 - 120 


0.3 X 10 8 


0.1 X 10~ 8 




0.7 - 0.8 


1.7 


CS 
















10 - 30 


10 x 10 6 


(0.001 - 0.3) X 


10 ~ 8 


0.2 - 0.5 


0.005 - 1.5 




35 - 40 


3 X 10 6 


0.01 X 10 -8 




0.2 - 0.3 


0.3 




45 - 80 


1 X 10 6 


0.1 X 10~ 8 




0.2 - 0.3 


0.05 




85 - 120 


0.3 X 10 8 


(0.3 - 1) X 10 


-8 


0.3 


0.009 - 0.03 


NGC6240: 














HCO + 
















10 - 25 


(1 - 3) X 10 4 


(0.3 - 1) X 10- 


8 


1.9 - 2.7 


0.4 - 0.7 




30 - 55 


1 X 10 4 


0.3 X 10~ 8 




1.8 - 1.9 


1.2 




60 - 120 


0.3 X 10 4 


1 X 10~ 8 




1.6 - 1.7 


0.7 


HON 
















10 - 20 


(0.1 - 1) X 10 6 


(0.3 - 1) X 10- 


7 


1.3 - 1.7 


0.03 - 0.1 




25 - 30 


0.3 X 10 6 


0.1 x 10~ 8 




1.3 - 1.4 


1.7 




35 - 55 


0.1 X 10 6 


1 X 10~ 8 




1.3 - 1.5 


0.3 




60 - 120 


0.3 X 10 8 


0.03 x 10~ 8 




1.4 - 1.5 


5.5 


t the chi-squarcd fit was 


calculated as — 




model) 


2 , where R a b s 


and R mo del are 



the observed and modeled line ratios, respectively, and is the associated error. 
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Table 8. Estimates of the total dense gas mass in Arp220 and NGC6240. See §4.41 for 

details. 





M denae \ 


X10 10 M Q ] 




Arp 220 


NGC6240 


Method 1 






HCO+(1-0) 


0.5-1.6 


1.3-2.2 


HCN(l-O) 


1.8-4.2 


1.0-2.8 


CS(2-1) 


0.3-1.5 




Method 2 






HCN 


1.4-7.0 


1.7-5.7 


Best Estimate 


0.3-1.6 


1.0-2.2 



